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ABSTRACT 

In Type la Supernovae (SNe la), the relative abundances of chemical elements are affected by the 
neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures 
in the spectra near maximum light, spectral features may be used to constrain the composition 
of the progenitor. We calculate the nucleosynthetic yields for three SN la simulations, assuming 
single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and 
calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and 
the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition 
scenario with different ^®Ni yields and the W7 simulation to control for differences between explosion 
models and total yields. While the overall yields of intermediate mass elements (16 < A < 40) differ 
between the three cases, trends in the yields are similar. With increasing metallicity, ^®Si yields 
remain nearly constant, '^°Ca yields decline, and Ti and ®^Fe yields increase. In the synthetic spectra, 
we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: 
a Ti feature around 4200 A and a Fe feature around 5200 A. In all three simulations, their pseudo 
equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two 
features may allow differentiation among progenitor metallicities of observed SNe la and potentially 
help reduce the intrinsic Hubble scatter. 

Keywords: nuclear reactions, nucleosynthesis, abundances - radiative transfer - stars: abundances - 
supernova: general 


1. INTRODUCTION 

Type la supernovae are generally considered to be the 
result of the thermonuclear disruption of carbon-oxygen 
white dwarfs in binary star systems. Although the path 
which leads to these explosions is uncertain, the results 
are easily seen and quite useful. PowCTed b y the radioac¬ 
tive decay of ®®Ni ( Colgate &: McKee|l969 ), SNe la have 
luminosities on the order of ~ 4 x 10'' Lq. This ex¬ 
treme luminosity in conjunction with a generally well 
follo wed correlation between peak brightn ess and decline 
rate (Phillips 1993 Phillips et al. 1999) allow SNe la 
to be used as standard candles out to redshifts of z«l. 
Measurement of the expansion history of the universe in 
this redshift range demonstrated that the expansion of 
the uni verse is accelerating due to the influence of dark 
energy (|Riess et al.||1998| |Perlmutter et al.||1999|). 

Though the 'I'ype la class as a whole is generally ho¬ 
mogeneous, variations betwee n individual events’ spec¬ 


tra and light curves do occur( Branch et al. 1993, 2006 
Blondin et al. 2012). Many of these variations can be 


traced back to the amount of iron-group material pro- 
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duced in the explosion with a portion of this material 
being radioactive ®®Ni. Explosions with lower amounts of 
iron-group material, and thus lower amounts of ®®Ni, are 
redder, dimmer, and have faster decline times than those 
with higher amounts of iron-group material. One factor 
that can affect the production of stable and radioactive 
iron-group material is the metallicity of the progenitor. 
Higher metallicity progenitors have larger abundances of 
neutrons causing greater amounts of stable ir on-group 
material to be produced and less radioactive Ni (Timmesj 

erM][2003t . 


Studies of explosion models used in this work (de¬ 
scribed below) confirm this result holds when considering 
the effect of the met allicity on the flame speed in addi- 


tion to composition (Townsley et al.|2009 Jackson et al. 


2010). The capability to observationally constrain the 


metallicity of SN la progenitors would reduce some of 
the uncertainty caused by intrinsic variations between 
individual events. In this study, we seek possible spec¬ 
tral features that could serve as indicators of progenitor 
metallicity. 

Observationally, it has been clear for some time that 
environm ent plays a ro l e in t he observed behavior of 
SNe la. Hamuy et al.| (2000) saw in a sample of 62 
SNe la that brighter, slower declining events preferen¬ 
tially occurred in spiral galaxies, while elliptical galax¬ 
ies were host primarily to quickly declining, dim events. 
Also, they find that metal-poor en vironments produced 
brighter SNe la. Kelly et al. (2010) find that SNe la that 
occur in larger, more massive galaxies were 10% brighter 
than other SNe la with similar light curve shapes. They 
claim that this trend could be caused by a correlation be¬ 
tween galaxy mass and metallicity, older progenitors in 
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higher mass galaxies, or o ther factors such as dust. Addi¬ 


tionally, Ellis et al. 
redshift SJNe la t 


(20081 found that for 36 intermediate- 
le ditierences between observed UV 


spectra could not be accounted for by dust alone requir¬ 
ing some other environmental factors such as metallicity 


to play a role. And, Howell et al. (2009a) found that 
SNe la that occured in host galaxies with higher metallic- 
ities produced up to 10% less ®®Ni. These observational 
studies all use host galaxy metallicity as a proxy for the 
progenitor metallicity. However as metallicity can vary 
by location and time in galaxies, it would be more effec¬ 
tive to have an indicator of the individual progenitor’s 
metallicity in the SN la spectrum itself. 


bration of SN absolute brightness and thereby decreasing 
the scatter from the Hubble law. In this present study, 
we do not attempt to resolve this broader question, but 
instead focus on evaluating, from theory, candidate in¬ 
dicators of the metallicity of the progenitor that are not 
too blended or obscured by neighboring features. Since 
metallicity is thought to be mostly a secondary param¬ 
eter relating to the SN brightness, the next step would 
be to attempt to separate it from the primary param¬ 
eter. Given our limited sampling of theoretical models 
and their current uncertainty, we refrain from addressing 
primary variation in brightness directly here, and instead 
limit our study to searching for candidate spectral indi¬ 
cators of metallicity. 

Past theoretical studies have focussed on metallicity ef¬ 
fects on synthetic spectra, in particular the UV flux, by 
making sensible ad hoc modific ations to the abund ance 


profiles of explosion models (see Brown et al. (2015), and 
references therein). Generally, such modifications are not 
necessarily self-consistent with different levels of progen¬ 
itor metallicity. For example, this type of modihcation 
does not account for the effect of progenitor metallic¬ 
ity on the ®®Ni production, which affects post explosion 
ejecta temperatures, and therefore the shape an d color 


of synthetic light curves and spectra. De et al. (2014) 
have shown that the production of intermediate inass el- 
ements is also affected by progenitor metallicity in a way 
that is expected to be fairly model independent. These 
changes in the yields of the intermediate mass elements 
(IMEs, with atomic mass number 16 < A < 40) should 
also have visible effects on the observed spectra. They 
find the greatest change to be in the abundance of Ga 
with the least amount of change in the abundance of Si. 
This indicates that the characteristic Sill line at 6150 A 
may remain fairly static with metallicity, but the fea¬ 
tures from other IMEs should vary in such a way that 
the progenitor metallicity could possibly be measured. 

Here, we use multi-dimensional simulations of SNe la 
using the deffagration-to-detonation transition paradigm 
with the FLASH hydrodynamics code in conjunction 
with particle post-processing and the 1-D radiative trans¬ 
fer code PHOENIX to examine the effects of progenitor 
metallicity on nucleosynthetic yields, light curves, and 
spectra. In section 2, we outline the SN la models we 
are using as well as a description of the hydrodynamic, 
particle post-processing, and radiative transfer calcula¬ 


lationship between spectra and brightness 

(e.g. iNugent 

variations of the 

et al. 1119951 IBailey et al.||20091 |Blondin et 

al. |2Ulla|bl 

simulation of a C 


tions. In section 3, we discuss the nucleosynthesis re¬ 
sults, synthetic lig ht curves, and spectra. We also com¬ 
pare our resu lts to Foley & Kirshner| (|2013 ) and Graham 
et al. (2015), who examined SIN2011 by and sKl2Ullfe 
two spectroscopically normal SNe la which had nearly 
identical decline times and optical spectra, but interest¬ 
ing differences in the NUV region. Finally, in section 
4, we give our conclusions and discuss the possibility 
of using spectral features as an indicator of progenitor 
metallicity. 

2. SN la EXPLOSION SIMULATIONS 

We consider two types of explosion simulations, both 
lelayed-detonation model for SNe la 
The main model is a 2-dimensional 
lhandrasekhar-mass WD exploding via 
the deflagration-detonation transition (DDT) mecha¬ 
nism. These simulations are identical to those we have 


this scenario ( 

Krueger et al. 

2010 

Jackson et al. 

2010 

Krueger et al. 

2U12I. tor comparison, we also compute 


well-studied, 1-dimensional W7 model (Nomoto et al. 


1984). We briefly describe our 2D explosion model and 


post-processing here. 

2.1. Hydrodynamic Simulations 
We chose two realizations of a GO white dwarf with 


a central density 2 x 10® g/cm®. Krueger et al. (2012) 
created 5 CO white dwarf progenitors with central den- 
sities that ranged from 1 x 10® to 5 x 10® g/cm®. For 
each progenitor, 30 realizations were created by apply¬ 
ing initial conditions in the form of random perturba¬ 
tions on the initial burned region. This was done to bet¬ 
ter understand potential systematic biases such as how 
the morphology of the initial conditions might influence 
the results. Many of their progenitors produced higher 
yields of ®®Ni than is typical in SNe la. In light of this, 
we selected realizations ROl and RIO due to their lower 
estimated ®®Ni yields, 0.8 Mq and 0.7 Mfr^ respectively, 


which are closer to typical values for SN la (Howell et al. 
2009b|). In this paper we refer to these two simulations 
as UDT-high and DDT-low. 

We recalculated the simulations because the origi¬ 
nal simulations did not feature tracer particles that are 
needed for this work (see below). For this work we used 
a newer version of the FLASH code, version 4.0. The 
simulati on results are nearly i dentical to the those pub¬ 
lished in Krueger et al. (2012), with some small changes 
due to changes to the hydrodynamics method in the in¬ 
tervening FLASH releases. The overall software consists 
of the publically available version of the FLASH software 
instrumenlj^ along with added component^ for efficient 
modeling m carbon-oxygen flames and detonations de- 
veloped during our previo us work d Galder et "^ 2007[ 


Townsley et al.||2(307 |2009] [Jackson et al.||2010[ |Krueger| 

et al.||2012p ' as well as soft ware for computing nu cleo- 

synthesis in post-processing (Townsley et al. 2015 P. In 
addition to the hydrodynamics itself, the current public 
FLASH software release includes the advection-diffusion 

® available from http://flash.uchicago.edu 
® available from http://astronomy.ua.edu/townsley/code 
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treatment used to model the propagation of the defla¬ 
gration. The FLASH modules implementing the carbon- 
oxygen burning modeP will also be included in a future 
public release of the FLASH software instrument. 

Our simulations begin at the end of the core convection 
phase of the runaway just after the ignition of a pro pa- 
gating thermonuclear flame (see Nonaka et al. 2012 for 
a discussion of this ignition process). We make the as¬ 
sumption that the flame is ignited in such a way that it 
spreads in the center of the star before plumes rise to 
lower densities. This is an analog of the “multi-point” 
ignition scenarios explo red by other auth ors in 3D (e.g. 


Seitenzahl et al. 


1 2012) suggests t 


2013 for an example). Nonaka et al. 
lat multi-point ignitions are likely dis¬ 
favored by nature. Nevertheless, as others have shown 
previously and we demonstrate again in this work, this 
does lead to an explosion that reproduces spectral obser¬ 
vations of normal SNe la fairly well. Our simulations fol¬ 
low the nuclear and hydrodynamic evolution of the flame 
within the star until it reaches a density of 10^'^ g cm“^. 


by a DDT mechanism ( 

Khokhlov|1991 

Poludnenko et al. 

2011 
of a 
there 

1. 'i’his detonation is ignited witi 
small heated region since all bn 
fore any physical DDT processes. 

1 the introduction 
rning fronts, and 
are unresolved in 


our full-star simulations. The complete incineration of 
the star by the detonation then ensues. Active burning 
fronts cease due to the expansion of the star by about 1.8 
seconds and then we follow the expansion of the star to 
4 seconds, at which point it is already close to force-free 
expansion. 

Our model of carbon-oxygen burning used in hydrody¬ 
namics is a three-level burning model that is calibrated 
against steady state deflagrations and detonations. The 
resulting temperature, T, and density, p, in the explosion 
simulation enable post-processing of Lagrangian tracer 
particles (see next section) to obtain abundan ces behind 


the reaction front with an accuracy of 5-10% (Townsley 
et al.|[20l5). The global energy deposition based on the 


burning model is consistent with the yields obtained from 
post-processing to within 10%. 

2.2. Lagrangian Tracer Particles 

During the explosion simulations, the temperature- 
density histories of a set of 100,000 Lagrangian tracer 
particles are recorded. These tracer particles are inac¬ 
tive, meaning that they do not affect the hydrodynamics 
but are only advected with the flow, recording the his¬ 
tory of local fluid variables. During the burning phase of 
the explosion, the histories are recorded every time step, 
giving the full time reso lution data which is ne eded for ac¬ 
curate post-processing ( Townsley et al.|20T5 |. The equal 
mass particles are distributed randomly. After the hydro- 
dynamic simulation is complete, we use the TORClj^ 
software instrument to compute the evolution of a set 
of 225 nuclides subject to the conditions recorded for a 
given fluid element. This post-processing step allows the 
computation of detailed nucleosynthetic yields without 
having to compute the reaction or advection of all of 
these species in the main hydrodynamic simulation. 

The density and temperature histories for W7 were 

Available at http://cococubed.asu.edu Our version available 
at http://astronomy.ua.edu/townsley/code 


obtained from Friedel Thielemann (priv comm). These 
consist of about 500 time steps for 175 zones out to 4.15 
seconds. These histories were treated in the same manner 
as Lagrangian trac ks from our 2D simulations. Similar 
yields were used by|De et al.|(|2014D to investigate the de¬ 
pendence of Si-group yields on metallicity in comparison 
to trends expected from quasi-equilibrium calculations. 

2.3. Nucleosynthesis 

The nucleosynthetic yields from the explosion simu¬ 
lations are determined by post-processing the histories 
of Lagrangian tracer pa rticles. The reaction network 


calculation, described in Townsley et al. (2015), deter- 
mines mass fractions for each particle. As described in 
Townsley et al. (2015), for fluid that is processed by the 
deflagration we reconstruct, in post-processing, a por¬ 
tion of the Lagrangian history in order to obtain a more 
realistic p{t), T{t) near the artificially thickened flame 
front. We do not do so for detonations or for the W7 
element histories, instead using the p{t), T{t) history di¬ 
rectly. After post-processing, radioactive products with 
half-lives smaller than 5 days, according to beta decay 
rates in TORCH, have their calculated abundances con¬ 
verted into their decay products. Three important decay 
chains of three isotopes are calculated during the radi¬ 
ation transport simulations: ^®Ni, ®^Fe, and ^®Cr. The 
decay chain of ®®Ni to ®®Fe is well understood to be crit¬ 
ical in the powering of SNe la light curves. ®®Ni has a 
half-life of 6.075 dayand transitions by beta decay to 
®®Co. ®®Co has a halhlife of 77.236 dayi^^ and transi¬ 
tions by beta decay to ^®Fe. The decay of ®^Fe to ®^Cr is 
complex. ®^Fe has a half-life of 8.275 hourjH^ and tran¬ 
sitions by beta decay to the 2-1- excited state of ®^Mn. 
The excited state has a half-life of 21.1 minuteJ^^^, but 
has two possible decay branches. 98.25% will transition 
by beta decay to ^^Co, and the remaining 1.75% decays 
to the ground state of ®^Mn via internal tr ansition. The 
ground state has a half -life of 5.591 dayJ^^ (jPessart et al.f 


2014 


Huo et al. 2007). However, we do not track this 


small amount of ®^Mn due to its negligible influence on 
spectral features. "‘^Cr has a half-life of 21.56 hourjHl 
and decays to The has a much longer half-life 
at 15.9735 day^, therefore its presence and changing 
abundance is important throughout the radiation trans¬ 
port calculations. 

The particles are then mapped onto a ID spherically 
symmetric velocity grid with 100 cells to match the map¬ 
ping of the radiation transport calculations, using the 
final velocities each particle has at the end of the sim¬ 
ulation. The mass fractions in each cell of the velocity 
grid follow from the average mass fractions of all parti¬ 
cles in that cell. These mass fractions per cell, together 
with a spherically symmetrized density profile from the 
hydro simulatiorp^ and the geometrical volume of each 
cell, give the totmyields (in gram) per species. 

If a velocity grid cell contains many particles we limit 
the number of particles that are post-processed to 100, 
selected randomly. This reduces the amount of time 
needed for post-processing. An independently selected 

http://www.nndc.bnl.gov/chart/ 

'i'his is done with Nathan Hearn’s QuickFlash analysis tools, 
using the latest checkpoint from the hydro simulation for which 
the ejected material has not begun to run off the grid. 
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set leads to abundances that differ on the order of 10% 
in velocity grid cells that contain greater than 100 parti¬ 
cles, and total yields that differ by less than 5%. 

2.4. Progenitor White Dwarf Abundances 

The progenitor white dwarf is comprised broadly of 
two regions: a convective core, the inner 0.8 Mq of the 
star, and a non-convective outer region. Our progeni¬ 
tor abundances are constructed to be similar to those in 


Dominguez et al. (2001), as well as to accoun t for pre- 
explosion carbon burning ( Piro fc Chang| 2008). 

In this work, the progenitor metallicity is varied by 
changing the initial abundance state of the tracer par¬ 
ticles, not by calculating separate hydro simulations for 
different initial abundances, because the energy release 
by carbon-oxygen burning is only marginally modified 
by the metallicity variations. The hydrodynamic initial 
model is at a single metallicity of Z/Zq = 1.33. The 
three parameters that control the carbon-oxygen burn¬ 
ing model used in the hydro stand for the abundances of 
and ^^Ne. We choose the convective core to be 
made up of 40% 3% ^^Ne, with the remaining 57% 

as the envelope is 50% 48% and 2% ^^Ne. 

The ^^Ne is a stand-in for the neutron excess contributed 
fro m the metallicity in b oth the core and the outer lay¬ 
ers (Timmes et al. 2003). In the convective region the 
additional 1% of ^^Ne stands in for the neutron excess 


from the products of pre-detonation simmering (Piro & 
Bildsten||2008 Chamulak et al.||2008). 

When post-processing the particles with a large nuclear 
network, more complete and realistic progenitor abun¬ 
dances can be used. Figure shows how the material 
affected by metallicity and the simmering ashes are dis¬ 
tributed through the star for both the reduced abundance 
set and the abundances used for post-processing. A par¬ 
ticle’s initial location is used to determine whether or not 
the particle began in the convection region, and thus the 
initial fraction and the relative contributions from 
metallicity and simmering. The initial abundances for 
post-processing are then set using 4 components: 
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Figure 1. Top: The initial composition and temperature (black) 
profile of the progenitor white dwarf used in the hydrodynamic 
calculations. This reduced set of nuclides consists of (red), 
(blue), ^^Ne (green). ^^Ne is a stand in for the contribu¬ 
tions from the simmering ashes (dashed), only present in the con¬ 
vective region, and the metallicity (dash-dot) present throughout 
the star. Bottom: Initial composition profile used in the post¬ 
processing showing a small subset of 225 species used, with values 
shown for different metallicities. (orange), ^®Ne (brown), and 
^^Ne (cyan) are all components of the simmering ashes held at 
equal fixed abundances in the convective region, and consequently, 
overlap one another. Outside of the convective region, ^®Ne also 
contributes as a component of the metals. 


The amount of simmering ashes, Xsim, is set to obtain 
consistency with the progenitor metallicity of ZjZ^ = 
1.33 used in the hydrodynamics, based on the Yg, 


i 


A,’ 


( 1 ) 


• The mass fraction Xi 2 q is set to the value at the 
particle’s initial position in the hydrodynamic pro¬ 
genitor. 


• The abundance of all other nuclides are first set to 
scaled solar abundances based on Z and abun¬ 
dances from Anders & Grevesse (1989). In doing 
this, we assume all metals lighter than ^®0, i.e. 

and are converted into ^^Ne and are 
thus added to that abundance instead. 


• A simmering contribution is added to some abun¬ 
dances for particles that are initially in the con¬ 
vective region. We approximate that the simmer¬ 
ing ashes, of a fraction Xgim determined below, are 
made up of an even mix of ^°Ne, and 

^^Ne ( Chamulak et al.||2008 ). 


• Once the ^°Ne, ^^Ne, ^^Ne and all heav¬ 

ier metal abundances are set in this fashion, the 
remainder of the material is assumed to be ^®0. 


For this study, the contribution from simmering ashes 
is held fixed as the metallicity is varied. Changing the 
metallicity thus affects the second and fourth step. 


in the outer layers. Here Zi and Ai are the number of 
protons and number of nucleons in species i. Matching Yg 
between the progenitor and the initial abundances used 
for post-processing gives 


^'^^Ojprog T A"i2(;; p,.Qg 


- V 2 


2Ne,prog 
6 


10 


= X16Q -h Vi2r;-h XzYe z 

16 12 ’ 


■ w 


1 

'4 V 16 


6 

I3 


10 10 


( 2 ) 


where abundances with the “prog” subscript are those 
found at the initial particle position in the progeni¬ 
tor used in the hydrodynamics, Xz is the sum of the 
scaled solar abundances of all material in the metallicity- 
scaled component, including ^^Ne, Yg.z is the Yg of the 
metallicity-scaled component, and equal abundances of 
the 4 simmering products have been assumed. For the 
hydrodynamic abunda nces given above, ZlZp) = 1.33, 
solar abundances from Anders & Grevesse (19^), and 
X 160 = 1 — Vsini — Xz — A 12 c, Fquationl^ gives 
Vsiin = 0.0175 in the convection zone and ATsim = 0 out- 
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side it. This value of Xsim is then used to construct the 
initial abundances for post-processing at all metallicities. 

2.5. Radiation Transport Simulations 
We calculate light curves and spectra using Phoenix. 


software instrument 1 

Hauschildt & Baron 

1999 

2004 

van Rossum 

2012 a 

). 

Phoenix numencail 

y solves the 


efficient and highly accurate short characteristic and op¬ 
erator splitting methods. It samples millions of atomic 
lines individually. Phoenix solves for the time evolu¬ 
tion using the radiation energy balance method. The 
code does not use the Sobolev approximation, diffusion 
approximations or opacity binning approximations. We 
operate Phoenix in one dimension, assuming spherical 
symmetry. 

The deterministic radiation transport method em¬ 
ployed in Phoenix calculates spectra without the noise 
that is inherent to Monte Carlo methods. This allows 
the study of small effects of individual atomic transi¬ 
tion line s on the calculated spectrum through knockout 
spectra (van Rossum|201^). Knock-out spectra are cal¬ 


culated by post-processing a normal calculation of light 
curves and spectra. For the calculation of knockout spec¬ 
tra, the opacity of individual transition lines is artificially 
set to zero while keeping the temperatures, electron den¬ 
sities, and atomic occupation numbers fixed to the values 
determined in the original calculation. The difference be¬ 
tween the original spectrum and a knockout spectrum is 
attributed to the same line opacity that was artificially 
set to zero. Here, as a diagnostic, we use knockout spec¬ 
tra to identify which chemical elements are responsible 
for the differences that metallicity changes cause in the 
calculated spectra. 

We note that a spectrum is not merely a linear combi¬ 
nation of individual lines. When multiple knockout spec¬ 
tra are shown together in one plot each of the knockout 
spectra is interpreted as described above, but there is no 
useful interpretation to a linear combination of knockout 
spectra. 


3. RESULTS 

3.1. Nucleosynthetic Yields 

The total yields for six different progenitor metallici¬ 
ties, Z/Zq e [0.1, 0.5, 1.33, 2.0, 3.0, 4.0], of each of the 
three explosion simulations are listed in Table 
The yields are grouped into Carbon and Oxygen (CO, 
A = 12 and A = 16), intermediate mass elements (IME, 
16 < H < 40), and iron-group elements (ICE, A > 40). 
56Ni is also included in the table because of its impor¬ 
tance as energy source for light curves and spectra. Ta¬ 
ble d] shows that the amount of radioactive Nickel de¬ 
creases with increasing progenitor metallicity while the 
total yields of iron-group elements increases due to the 
increased neutronization. 

The relative change in abundances with progenitor 
metallicity is plotted in Eigure The effect of pro¬ 
genitor metallicity on the nucleosynthetic yields of the 
explosion are fairly similar among the three models for 
most elements, especially for the DDT-high and DDT- 
low models. The most marked difference between W7 
and the 2D DDT models is that the Ne and Na yields 
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Figure 2. Ratios of the nucleosynthetic yields of high (Z=2 Z©) 
versus low (Z=0.5 Z©) progenitor metallicities for the DDT-high 
model (top), the DDT-low model (middle), and W7 (bottom) on 
a logarithmic scale. A value close to 0, like Si, means that the 
same nucleosynthetic yields were obtained for different progeni¬ 
tor metallicities. The effect of progenitor metallicity on the total 
nucleosynthetic yields of the explosion are fairly similar between 
the two DDT models. The affects on W7 are fairly consistent in 
sign with the DDT models except for a number of elements above 
scandium. 















































































































































6 


Table 1 

Nucleosynthetic yields grouped by atomic mass in Mq 
for three simulations and six progenitor metallicities. 


Z/Zq 

0.1 

0.5 

1.33 

2 

3 

4 


DDT-high 

o 

o 

0.053 

0.053 

0.053 

0.053 

0.052 

0.052 

IME*^ 

0.322 

0.320 

0.314 

0.310 

0.302 

0.295 

IGE'= 

0.992 

0.994 

1.003 

1.006 

1.012 

1.025 


0.848 

0.831 

0.800 

0.775 

0.738 

0.703 


DDT-low 

CO 

0.076 

0.077 

0.077 

0.076 

0.076 

0.075 

IME 

0.385 

0.382 

0.375 

0.369 

0.360 

0.351 

IGE 

0.906 

0.909 

0.916 

0.922 

0.932 

0.941 


0.771 

0.756 

0.726 

0.703 

0.669 

0.637 


W7 

“CD 

0.166 

0.166 

0.165 

0.163 

0.161 

0.159 
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Figure 3. Total yields in Mq of ^®Si (red), ®^Fe (green), 
(black) as a function of metallicity for the two realizations, DDT- 
high (circles, solid line) and DDT-low(squares, dashed lines), as 
well as W7 (triangles,dash dot lines) 




Figure 4. Ratio of ®'^Fe to ^*Si yields (top) and '*‘^Ca to ^®Si yields 
(bottom) from the 2D DDT and the W7 simulations. The ratios 
of these abundances show similar smooth trends with progenitor 
metallicity for these three models. 


increase more with metallicity in the W7 model. 

Figure|^shows the integrated yields of ^®Si, '^°Ca, and 
^'^Fe for the six progenitor metallicities. We find that 
with increasing metallicity, the overall yields of ^®Si stay 
relatively constant. Conversely, the total yields of ^°Ca 
decrease and ®‘^Fe increase with increasing m etallicity 
similar to the trends seen in De et al. (2014) with de¬ 
creasing Ye. The relatively flat trend in the ^®Si yields 
makes this species a useful reference for comparing the 
trends in the yields of other species between the three 
models. The ratios of '^^Ca and ®‘*Fe relative to ^®Si are 


shown in Figure]^ As expected from Figure!^ the trends 
are similar among the models, though the absolute ratio 
of ^°Ca to ^®Si is different in the W7 model due to Ca 
being produced in higher amounts at a wider extent in 
the DDT models (Figure]^. The trends are also smooth 
over the metallicity range that we are studying. 

Figure shows the chemical abundance profiles in ve¬ 
locity space from the DDT-high, DDT-low, and W7 mod¬ 
els for two metallicities: Z=0.5 and 2.0 Zq. Differences 
in explosion model produce differences in the profiles for 
a single metallicity. While all three models feature a ®®Ni 
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Figure 5. Chemical abundance profiles in ID spherical velocity space of the DDT-high model (top row), the DDT-low model (middle 
row), and the W7 model (bottom row) for two different initial metallicities, Z=0.5 and Z=2. 



















































core that extends to approximately 1.0 x 10® cm/s, W7 
is the only model that has a “nickel hole”, which is a 
known characteristic of 1-D models. In the higher ve¬ 
locity portion of the ejecta, from « 1.6 x 10® cm/s to 
2.5 X 10® cm/s, the W7 model’s ejecta consist of unburnt 
C, O, and Ne, much simpler than that produced in the 
DDT-high and DDT-low models. In this region of the 
material, the DDT-high and DDT-low models’ profiles 
consist of predominantly C and O, with < 10% of Si, 
Mg, Ne, S, Ni, Ar, and Ca. At the highest velocities, 
> 2.7 X 10® cm s“^, the DDT-high and DDT-low ejecta 
take on a similar configuration to that of the outer layers 
of W7, unburnt C, O, and Ne. At intermediate veloci¬ 
ties, between 1.0 x 10® cm/s and 1.6 x 10®, Si peaks in 
all three models. This is also the location of other IME 
peaks such as S, Ar, and Ca. 

Examining the same explosion models at different 
metallicites presents additional features. At Z = 2 Zq, 
®®Ni makes up less of the total Ni yield in the three ex¬ 
plosion models, though there is little change in the over¬ 
all mass fraction of Ni. Overall, the Ee mass fraction 
increases with the increased metallicity. However, not 
only is the mass fraction of Ee higher, the Ee rich layers 
extend further outward into the higher velocity ejecta. 
In all three modes, the Ee rich layer reaches into the 
material that is mostly intermediate mass elements. The 
increase in metallicity also causes the Ca peak to be both 
lower and narrower. The Si peak is largely unaffected by 
the increase in metallicity. This is also reflected in how 
little the overall Si abundance varies in Figure and the 
small change in the Si ii spectral feature at 6150 A with 


metallcity as shown in section 3.3 
Figure shows the abundance proHles of Si and Ti 
in the DDT and W7 models 30 days post explosion. 
It demonstrates how the Ti production is sensitive to 
the progenitor metallicity while the Si production is 
largely unaffected. The variation of Ti with metallicity is 
strongest in the region of incomplete Si burning and the 
unburned regions. A shift in the quasi-equilibrium abun- 
danc es due to metall icity is expected in the Si-burning re¬ 
gion ( De et al.|2014 ). The direct Ti production is mainly 
in the form of '’‘'‘Ti, which has a half-life of 60 yearj^ 
while the progenitor metallicity contributes mainly ^®Ti. 
Any '^®Cr produced during the explosion will decay with 
a half-life of 21.56 hour^Slfo which will decay fur¬ 
ther with a half-life of 15.97 dayJ^ to stable "^^Ti. We 
show the profile at day 30 post explosion to facilitate 
disc ussio n of Ti spectral features at this epoch in Sec- 


tion|^ At this time, 71% of the ^®Cr produced during 
the explosion has decayed to ^®Ti so that the Ti curves 
in Figure in some regions are higher than the “direct 
Ti” curves. 


3.2. Light Curves 

The bolometric light curves for the three explosion 
models with two extreme metallicities are shown in Fig¬ 
ure Higher progenitor metallicities systematically 
make the bolometric light curve rise more slowly and 
peak lower. The slower rise is caused by the higher 
abundance of metals in the outer layers of the ejecta, 
which raise the opacity, and make radiation diffuse out 



DDT-low 



Velocity [cm/s] le9 


W7 



Velocity [cm/s] le9 

Figure 6. Radial abundance profiles of Si and Ti at day 30 post 
explosion in the DDT-high model (top), DDT-low model (middle), 
and the W7 model (bottom). Different line colors represent dif¬ 
ferent initial metallicities. The Ti yield (solid line) in the outer 
regions of the burning zone is sensitive to the progenitor metal¬ 
licity, while the Si production (dashed line) is largely unaffected. 
Also shown, for each metallicity, is the portion of the Ti that is not 
a product of ^*Cr decay (direct Ti, dotted lines). 


http://www.nndc.bnl.gov/chart/ 
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DDT-high 




Figure 7. Bolometric light curves of the DDT-high (solid lines), 
DDT-low (dashed lines), and W7 (dotted lines) explosion simula¬ 
tions at two extreme progenitor metallicities of 0.1 (blue lines) and 
4 times solar (green lines). The light curves of explosion simula¬ 
tions with high metallicity progenitor are brighter at peak than 
those with low metallicity progenitors because the former have a 
higher ®®Ni yield. 


more slowly from the hot ®®Ni-rich core. The lower peak 
brightness is caused by the lower total mass in the 
models with higher progenitor metallicities, as shown in 
Table HI 

In Figure shows the UBVRI-band model light 
curvescompar M to the light cur ves from the data driven 


DDT-low 


model MFCS Jha et al. (2007). The systematic trend 
with progenitor metallicity in most light curves is small 
but significant, especially in the UBV-bands. This is 
apparent in Figure where the light curve shape pa¬ 
rameters, peak brightn ess and decline rate A Mi^fB), are 
plot ted, along with thelP hillips et al. (1999) relation (see 
also Mandel et al. 201^ i'he mod^ light curves have 
slightly brighter peak niagnitudes than the bulk of ob¬ 
served SNe la at the same /S.mi^{B). We find for the 
two 2D DDT cases used here, which have above-average 
yields, that the changes in ignition conditions that lead 
to different ®®Ni yields caus e the shape paramet ers to 
change perpendicular to the Phillips et al. ( 1999) rela¬ 
tion. This c onfirms a similar result found by peitenzahl 
et al. (2013) using 3D DDT models. This failure of mul¬ 
tidimensional DDT models to reproduce the Phillips re¬ 
lation in the expected way makes it difficult to make 
comparisons with empirical metallicity-luminosity rela¬ 
tions. In contrast, the progenitor metallicity introduces 
a change in the B-band light curve shapes that roughly 
follows the observed width-luminosity relation. 

The bottom panel of Figure shows the bolometric 
peak magnitude and decline rate. While higher metal¬ 
licity leads to lower bolometric peak magnitudes, just as 
it leads to lower S-band peak, the effect on AMi 5 (bol) 
is very small. This indicates that the range of AMi5(i3) 
(top panel of Figure spanned by metallicity is due 
entirely to a color eff^, in which the B-band light de¬ 
creases more quickly at higher metallicity while the de¬ 
cline rate of the total energy output is mostly unchanged. 
AMi 5 (bol) is even slightly slower at higher metallic¬ 
ity, the opposite of AMi^IB). In reality, any color ef- 




Figure 8. Bolometric and UBVRI-band light curves of the DDT- 
high (top), DDT-low (middle) and the W7 model (bottom), com¬ 
paring progenitor metallicities of Z=0.5 (solid lines) and Z=2 
(dashed lines). As a reference, MLCS light curves are plotted for 
a value of MLCS A = 0 (dash-dotted lines). The models are simi¬ 
lar in brightness and color evolution, especially the DDT-high and 
DDT-low simulations. The biggest difference is the peak brightness 
in the U-band between the 2D DDT simulations and W7. 
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Figure 9. Peak B-band magnitude (Mb) vs. B-band decline rate 
(AMi 5 (i?)) {top panel) and peak bolometric magnitude vs. bolo- 
metric decline rate {bottom panel), for the DDT-high, DDT-low, 
and W7 models with initial metalliciti es of 0.1, 0.5, 1.33, 2, and 4 
times solar. The |Phillips et al.| ||1999|l relation is plotted in gray, 
with the shaded regions indicating the 1 and 2 cr observed variation 
also from that work. The range in AM 15 is present in the B-band 
but not the bolometric light curves showing tha t this is a color 
effect, as pointed out by |Kasen fc Woosley| p007|l 


fects caused by progenitor metallicity would be convolved 
with environmental factors such as extinction from dust. 



progenitor metallicity. 


3.3. Spectra 

Model spectra, calculated for W7 and the two DDT 
simulations, with different progenitor metallicities, and 
at fo ur d ifferent epo chs post explosion, are shown in Fig¬ 
ures through [ 1 ^ 

The mode l spectra are compared against the |Hsiao| 
et al. (2007) spectral templates. These templates rep¬ 
resent the typical spectral shape of normal SN la, as 
they are constructed as an average over many observed 
spectra. We set day 0 of the Hsiao07 templates (which 
is the time of maximum brightness) to the time that the 
B-band luminosity peaks in the model spectra, averaged 
over the set of progenitor metallicities, and fit a “dis¬ 
tance” to the imaginary source of the Hsiao07 templates 
to simultaneously fit the spectra of all epochs. We do 
not apply any other corrections, such as stretch or color 
corrections, to the Hsiao07 templates. 

The model spectra generally agree reasonably well with 
the templates, but there are a few significant discrepan¬ 
cies. First, none of the model spectra reproduce the far 
UV flux below 2500 A as it is observed on Day 10 and 20 


post explosion (pe). Second, although the Sill feature 
at 6150 A is present in the model spectra none of them 
accurately reproduces the observed shape of this feature. 
It is unclear whether this is a consequence of an under 
or over produced nucleosynthetic yield or a missing line 
or combination of lines from the radiative transfer calcu¬ 
lation. Third, the W7 model spectra do not reproduce 
the observed strong absorption feature around 3300 A at 
day 10 pe. 

Different progenitor metallicities give rise to variations 
in the model spectra. These variations change with time. 
Generally, the variations appear to affect the flux more 
on the blue side of the spectrum than on the red side 
and, over time, the range in which variations are strong 
gradually extends further to longer wavelengths. At day 
10 pe, the DDT and W7 spectra are most effected in the 
2600-3600 A range. At days 20 and 30 pe that window 
extends to 2600-4500 A and 2600-5200 A, respectively. 
In the latest spectra, at day 40 pe, the effect of progenitor 
metallicity presents only a small amount of variation. 

The systematic decrease of blue (short-wavelength) 
flux with progenitor metallicity is correlated in part to 
the decline in total ^®Ni yield with metallicity. This 
makes it more difficult to use this region as an indicator 
of metallicity that is independent of overall ^®Ni yield. 
The total ®®Ni yield affects the brightness of the light 
curve as well as the temperature of the radiating ejecta, 
and thus the spectral color. At the same time, the pro¬ 
genitor metallicity affects the composition of the SN la 
ejecta, which may affect the strength of individual spec¬ 
tral features. 

In order to help demonstrate the difference be¬ 
tween the ®®Ni-mass-temperature-color effect and the 
composition-feature-strength effect, Figure [l4| compares 
spectra at similar ®®Ni yields. From the top panel spec¬ 
tra at days 10 and 20 pe, we see that the region of the 
spectrum around 3000 A does not show as strong a metal¬ 


licity variation in this comparison as it does in Figure 10 
The case in the bottom panel with slightly higher yielc 
and lower metallicity, shows more metallicity variation 
in this region. Together these comparisons show that it 
will be necessary to carefully control for total ®®Ni yield 
in order to use a metallicity indicator from this spectral 

region. _ 

Features that vary consistently in both Figures 
and are good candidates for robust metallicity in¬ 
dicators that may not require careful control for ®®Ni 
yield, which renders the feature useless if the ®®Ni yield 
is not known. Two features that clearly vary with pro¬ 
genitor metallicity at day 30 pe in Figure [T3| as well as 
in all spectra in Figure 10 are 1) an absorption fea¬ 
ture around 4200 A, and 2) an absorption feature around 


5200 A. These two features are potential spectral indi¬ 
cators of progenitor metallicity. We will examine these 
features more below. 

There are a couple of other features that are promising, 
but their variation with metallicity is less con sistent for 
both varying (Figure [To| and fixed (Figure [l4|) ®®Ni. The 
variation with metallicity around 4400 A at day 20 pe in 
Figure is much less clear in Figure The feature 
around 5500 A at day 20 pe that appears to vary in Fig¬ 
ure does not vary in the same way in Figure |10| For 
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Figure 10. Model spectra of the DDT-high model (top) and DDT-low model (botto m) for five different progenitor metallicities at day 
10, 20, 30, and 40 post explosion (pe). For comparison against observational data, the |Hsiao et al.| | |2007[ | spectral templates are plotted 
in grey. The spectral feature in the black box around 4200 A and 5200 A at day 30 pe are potential spectral indicators for the progenitor 
metallicity (see text). 
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Figure 11. Spectra of Figure [To] extended to the IR range 
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Figure 12. Model spectra of the W7 mod el for five different p re-explosion metallicities at day 10, 20, 30, and 40 post explosion (pe). For 
comparison against observational data, the [Hsiao et al.| | |2007| | spectral templates are plotted in grey. The spectral features in the black 
boxes around 4200 A and 5200 A at day 30pe are potential spectral indicators for the progenitor metallicity (see text). 
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Figure 13. Spectra of Figureextended to the IR range 
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M(Ni56) = 0.70 Msun 



M(Ni56) = 0.77 Msun 

le40 



Wavelength [ang] 

Figure 14. Model spectra at similar ®®Ni yields, 0.7 Mq (top) and 0.77 Mq (bottom), from different DDT models at different progenitor 
metallicities. This comparison, together with Figure |10| helps to determine which variations in the spectra are due to changes in the 
abundance profiles of species other than ®®Ni. The most characteristic feature that robustly varies with progenitor metallicity, even at 
constant ®®Ni mass, are the absorption features highlighted with black boxes around 4200 A and 5200 A at day 30 pe. 
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Figure 15. The pseudo equivalent width (pEW) of the spectral 
features around 4000-4400 A (solid lines) and 4900-5400 A (dashed 
lines) as measured in the synthetic spectra of three explosion mod¬ 
els DDT-high (Mue). DDT-low (green), and W7 (red) at day 30pe 
(see Figures |10| and |14[ |. Using the three models investigated in 
the paper, these two absorption features gradually become stronger 
with increasing progenitor metallicity, making them good potential 
candidates for spectral indicators of SN la progenitor metallicity. 


this latter feature it appears that the two DDT mod¬ 
els produced different amounts of the chemical elements 
that are responsible for these features, and so this is not 
a simple effect of progenitor metallicity. 

The Sill P-Cygni feature at 6150A is present in the 
DDT and W7 model spectra at peak brightness and be¬ 
fore. Even though, as noted before, the observed shape 
of this feature is not very accurately reproduced in the 
model spectra, the change in metallicity has very little 
effect on both the depth of the absorption wing and the 
height of the emission wing. This is in agreement with 
the observation in Figures]^ and namely that the Si 
abundance remains relativmy static by changes in the 
progenitor metallicity. This allows for its use as a pos¬ 
sible calibration when examining the strengths of other 
features. 


3.4. Analysis of Spectral Features Around 4200 A and 

5200 A 


Two spectral features were identified above as potential 
spectral indicators of progenitor metallicity, in all three 
models: 1) an absorption feature around 4200 A, and 2) 
an absorption feature around 5200 A; both at day 30 pe. 
Both of these features become systematically stronger 
with increasing progenitor metallicity. Figure shows 
the trend of the feature strength with progenitor metal¬ 
licity, with the pseudo equivalent width pEW as measure 
for feature strength. 


pEW 



-FpCont(A)-E(A) 

F(A) 


dA, 


( 3 ) 


where FpCont is the pseudo continuum flux level as ap¬ 
proximated by a linear interpolation between the flux 
levels F at beginning and end points of integration Aa 
and A;,. The interpolation range for the two features is 
chosen as [4000, 4400] and [4900, 5400] in units of A. 


While the dependence of the pEW on metallicity 
is promising, the variation across metallicities in the 
ZjZQ =0.5 to 2 range is similar to that among the three 
explosion models considered her e. This means that, in 
order to use a plot like Figure[^to determine the metal¬ 
licity, it is necessary to obtain a control variable that dis¬ 
tinguishes among the models. Unfortunately the most 
accesible control variable Ami 5 (i?) does not appear to 
fill this role well. The DDT-high model at ZjZ^ = 0.5 
and the DDT-low model at ZjZ^ = 2 have both very 


similar Ami^{B) and very similar 4000-4400 A pEWs. 
Our limited sample of explosions available here does not 
allow us to robustly test control parameters. However, 
we are hopeful that a larger sample, theoretical or ob¬ 
servational, may be able to exploit the trends in pEW 
shown here to calibrate a measure of metallicity. 

In order to investigate which chemical species are re¬ 
sponsible for these two specific characteristic features we 


use knockout spectra (see Section 2.5). Figures 16 and 17 
displays a series of knockout spectra tor the DDT-fow am 
the W7 model at day 30 pe for Z/Zq = 0.1, 1.33 and 4. 


These knockout spectra are created by removing the 
opacity from all transition lines of a single chemical ele¬ 
ment from the radiative transfer calculation and repeat¬ 
ing the simulation with fixed temperatures, electron den¬ 
sities, and population numbers as in the original calcu¬ 
lation. The recalculated spectra show the effect of the 
absence of opacity. This allows us to identify what and 
how strongly different elements are involved in forming 
specific spectral features. 

In the knockout spectra. Figures 16 and [m the full 
spectra are represented by the solid mack line. The color 
patches represent the knockout effect that each element 
has on the full spectrum. Patches below the spectra indi¬ 
cate emission features: the flux level is that much lower 
in the absence of opacity from that species. In the same 
way, patches above the full spectrum represent absorp¬ 
tion features. 

Using the knockout spectra, we see that the two identi¬ 
fied spectral indicators of progenitor metallicity, around 


4200 A and 5200 A at day 30 pe, can be associated with 
Titanium and Iron production, respectively. The higher 
the progenitor metallicity the stronger the Titanium and 
Iron abosption features appear in the spectra. Iron also 
shows a shift from emission to absorption in this feature 
as metallicity increases. Looking back at Figures]^ and 
we see that indeed the Ti and Fe abundances increase 
with the progrenitor metallicity in the outer regions of 
the ejecta, i.e., at velocities > 10® cm/s. 


3.5. Spectra Ratios 


Foley fc K irshner (20131, and more recently, Graham 
et ai.|(20I5|) investigated the two SNe la SN 201 Ite and 
SN 2011 by. These two objects exhibited remarkably sim¬ 
ilar optical spectra and decline times, and even peak 
brightnesses after accounting for possible errors in the 
measured distance to 2011by. Thus the two objects be¬ 
come nearly identical. However, the two objects dis- 
play differences in the near UV regions of their spectra. 
Graham et al. (20151 propose that one possible source 
of these differences could be a difference in progenitor 
metallicity. To investigate this proposal, they first nor- 
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Figure 16. Knock-out spectra for the DDT-low model with metallicity Z/Zq =0.1 (top), 1.33 (middle), and 4.0 (bottom) at day 30pe. 
Contributions from line opacity of individual elements to the full spectrum (black solid lines) are represented by colored patches. These 
are calculated by removing the line opacities from a single element and recalculating the spectra while keeping the temperatures and 
population numbers fixed. Emission contributions are represented by colored patches below the full spectra, and absorption contributions 
are represented by patches above the full spectra. The spectral features highlighted with black boxes are two potential spectral indicators 
for progenitor metallicity. 
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Figure 17. Like Figure [T^ but for the W7 model. 
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malize the two spectra across 4000-5500 A and then take 
the ratio of the two spectra, llfe:llby. This ratio is 
then compared to ratios of model spectra a t difference 
metallicites produced by Lentz et al. (2000). By com¬ 
parison, they hnd the best match between the observed 
spectral ratio and the model spectra correspo nds to a 


ent blueshifting of features in the DDT-low spectra com¬ 
pared t o th e same features in the DDT-high spectra. In 
Section 13.31 


factor of 30 difference in progenitor metallicity (Graham 
et al.|20i5 Figure 8). In this section, we discuss how well 


our inodels reproduce the observed UV spectra in these 
cases, and challenges in making this type of comparison 
when metallicity and explosion strength (®®Ni mass) are 
allowed to vary independently. 

In Figure [TSl we take 4 pairs of spectra from our mod¬ 


els at day 20 pe a nd follow a similar procedure as Gra- 
ham et al.| (|2015|). Two pairs, DDT-high Z = 0.5 Zq 
& DDT-lowZr= 2.0 Zq and DDT-high Z = 2.0 Zq 
& DDT-low Z = 4.0 Zq, have similar decline times as 
shown in Figure The other two pairs, DDT-high Z 
= 2.0 Zq & DDT-low Z = 0.1 Zq and DDT-high Z = 
4.0 Zq & DDT-low Z = 4.0 Zq, have similar optical 
spectral features at day 20 pe. Each of the spectra are 
normalized across the 4000 - 5500 A region. The top 
and middle panels of the left column of Figure show 
the normalized spectra in log space of each pair selected 
for similar decline times. The top and middle panels of 
the right column of Figure [T8| show the normalized spec¬ 
tra in log space of each pair selected for similar optical 
features. Also included in the top and middle panels of 
both columns of Figure are the spectra of SN20 11fe 
at -2 days from peak brightness (Mazzali et al. |2014[ ) and 
SN2011by at -I day (both Melissa Graham, priv. coiiim.). 
These spectra were first binned in increments of 5 A and 
then normalized in the same way as our calculated spec¬ 
tra. Our calculated spectra share a similar sequence of 
features, especially blueward of about 3000 A, with these 
two observed spectra. However, the features in the model 
spectra are both stronger and located at slightly shorter 
wavelengths. Also the overall relative flux level blueward 
of about 2700 A is about a factor of 10 lower in the model 
spectra. 

The bottom panels of both columns in Figure [T8| show 
the ratios taken between each pair as well as the ratio of 
SN2011fe to SN2011by. One of the most prevalent differ¬ 
ences in the model spectra is in the 3000 A region, which 
is one of the regions most strongly effected by varying 
metallicity. The model pairs selected to have the same 
decline time show a factor of 4-6 ratio in this region, 
while the observed pair shows no significant difference. 
Gomparing the ratios of pairs, DDT-high Z = 0.5 Zq 
& DDT-low Z = 2.0 Zq (Figure 18 bottom left panel, 
red curve) and DDT-high Z = 4.0 Zq & DDT-low Z 
= 4.0 Zq (Figure 18 bottom right panel, green curve), 
demonstrates the difficulty in separating metallicity from 
®®Ni yield variations. The location of the peaks in these 
two pairwise comparisons is similar, having ratio peaks 
at about 1900, 2350, 2500, and 3600 A. However, while 
the first pair are at different metallicities, the second pair 
are at the same metallicity but with different ®®Ni yields. 
A possible factor in the similarity can be seen in the 
®®Ni yields. Looking at Table a it is clear that none of 
the four cases produced equal amounts of ®®Ni, though 
the ratios of ®®Ni produced in each pair are very similar 
at « 1.1. Another contribution comes from the appar- 


we have discussed that many of the changes 
in bluer side of the NUV region are caused by tempera¬ 
ture effects from the changing ^®Ni yield. Since changes 
in spectral features from this temperature effect occur 
alongside changes from the varying ejecta composition, 
isolating the two effects by taking spectral ratios appears 
challenging. A more conclusive comparison will require 
model spectra that are more similar to the observations 
in this spectral range. 

4. CONGLUSIONS AND DISGUSSION 

We have utilized a combination of multidimensional 
hydrodynamics calculations, nucleosynthetic particle 
post-processing calculations, and radiative transport cal¬ 
culations in a preliminary investigation into the role that 
progenitor metallicity plays in the chemical composition 
of SN la ejecta, and the consequences that has for their 
light curves and spectra. These calculations were per¬ 
formed on two 2-dimensional DDT models, as well as on 
the 1-dimensional W7 model. Lagrangian tracer particles 
were produced for the two DDT models using the FLASH 
hydrodynamics software instrument. These tracer par¬ 
ticles as well as the W7 temperature-density histories 
were then post-processed using TORGH with a 225 mem¬ 
ber nuclear reaction network to calculate nucleosynthetic 
yields and ejecta abundance profiles. We have used the 
PHOENIX radiation transport software to compute light 
curves and spectra for these ejecta profiles. With this ap¬ 
proach we examine the effects of progenitor metallicity 
on the nucleosynthetic yields of the explosion and on the 
synthetic light curves and spectra. 

Our approach is different from past theoretical stud¬ 
ies that have focussed on the effects of metal content in 


the ejecta on synthetic spectra (see Brown et al. (2015), 
and references therein) in two ways. I'irst, we calculate 
self-consistent nucleosynthetic yields from explosion sim¬ 
ulations by changing the progenitor metallicity instead of 
making modifications to post-explosion abundance pro¬ 
files. This takes into account any effects of progenitor 
metallicity on the ^®Ni production, which affects post ex¬ 
plosion ejecta temperatures, and therefore the shape and 
color of synthetic light curves and spectra, and in partic¬ 
ular in the UV. Second, instead of focussing on the UV 
flux, we explore the effect of progenitor metallicity on the 
strength of spectra features in the optical. The UV flux 
has been found to be sensitive not only to metal content 
in the ejecta, but also to other properties li ke the density 
struc ture and the type of explosion model (Brown et al.j 
2015 and references therein), a hnding that is conhrmed 
also in this paper. 

We Hnd that two aspects of SN la spectra present chal¬ 
lenges to finding spectral indicators of progenitor metal¬ 
licity or neutron excess. The first is that SN la spectra 
are formed by a large number of blended lines from a va¬ 
riety of species. This makes isolating individual features 
challenging, more so in some spectral ranges than oth¬ 
ers. The second challenge is that while the abundances of 
both the unburned material and the IME material vary 
systematically with metallicity, as expected, the result¬ 
ing variation in spectral features is confused with the 
variation of the total ®®Ni yield with metallicity. Our 
inclusion of explosions with different yields allows us to 
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Figure 18. Top and middle: The normalized maximum light NUV spectra of pairs of spectra chosen for their similar decline times (left) 
and similar optical spectral features at maximum brightness (right). The spectra of SN2011fe and SN2011by are also shown for comparison. 
Bottom: Ratios of the maximum light NUV spectra of each pair chosen for similar decline times (left) and similar optical spectral features 
(right) with the ratio of SN2011fe to SN2011by included for comparison. Ratios are shown for the top panel pair (red) and middle panel 
pair (green). 


explore this ambiguity. Hig her metallicity, whic h pro¬ 
vides a larger neutron excess ( Timmes et al.|2()03 ), leads 
to a lower yield of ®®Ni, instead favoring stable Fe-group 
yields. This reduction in the energy source for the pho- 
tospheric phase causes a lowering of the ejecta temper¬ 
ature, and reduces the overall flux in the UV. Isolating 
individual features, which vary due to changes in yields 
with metallicity, from this overall shift in flux intensity 
is challenging, especially in the early-time UV spectrum. 

Despite these challenges, we find two spectral features 
that show promise as potential robust spectral indica¬ 
tors of progenitor metallicity. These features stand out 
as metallicity indicators even when comparing explosions 
with similar total ®®Ni, a comparison for which many 
other features are the same (see Figure 14). The most 
promising feature we have found is a Ti feature at day 30 
post explosion that extends from 4000-4400 A, based on 
the knock-out spectrum. The pEW for this feature (Fig¬ 
ure 151 shows a robust and mostly consistent variation 


with metallicity for each of the three explosion simula¬ 
tions studied here. However, a control for overall yield 
will still be necessary for this to be used as a precision 
indicator of metallicity. The small sample of explosion 
models computed here does not allow evaluation of a 
suitable control variable, but a broader set of computed 
cases or an observational collection of spectra may. 

Our 2D DDT models of SNe la produce spectra that 


are sim ilar to those of observed SNe la in many respects 
(Figure [l0|), in some regions better than the spectrum of 
W7 (Figure [I^. It is important to note that we are not 
fitting spectra, so that a precision match is not the goal, 
rather similar overall spectral features. One of the main 
differences between the 2D DDT and W7 is apparent in 
Figurej^ where the DDT models have intermediate mass 
and iron-group material extending out to higher veloci¬ 
ties than W7. This appears to make the absorption and 
emission features in the 2D DDT spectrum, especially in 
the blue portion of the spectrum, more similar to those 
in observed SNe la. 

There are several trends in yields that result from 
changing progenitor metallicity. Table shows that for 
all three models, increasing metallicity causes a decline 
in the production of intermediate mass elements and an 
increase in the production of iron-group material. How¬ 
ever, as is expected with the increased neutronization, 
the increased production of these stable iron-group ele¬ 
ments comes at the cost of reduced production of ®®Ni. 
Also at higher metallicities, the iron-group material pro¬ 
duced in incomplete Si burning extends further out into 
the IME layers of the ejecta. The combination of the 
increased opacity from the iron-group material and the 
reduced production of ®®Ni results in light curves that 
are both slower rising and dimmer than those of lower 
metallicities. However, we find only an ~ .10 mag de- 
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pendence of the peak luminosity between the lowest and 
highest progenitor metallicities, and with our small sam¬ 
ple size we are hesita nt to make st r ong c omparisons with 
observations such as Kelly et al. (2010). Although the 
overall yield of interniediate mass elements declines with 
increasing metallicity, the yield and distribution of Si 
remain relatively static. Consequently, the SN la char¬ 
acteristic Si P-Cygni feature at 6150 A remains mostly 
unchanged with metallicity. 
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